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ABSTRACT 

Isophotal analysis of M87, using data from the Advanced Camera for Surveys, reveals a projected 
displacement of 6.8 ± 0.8 pc (^ O'.'l) between the nuclear point source (presumed to be the location of 
the supermassive black hole, SMBH) and the photo-center of the galaxy. The displacement is along a 
position angle of 307 ± 17° and is consistent with the jet axis. This suggests the active SMBH in M87 
does not currently reside at the galaxy center of mass, but is displaced in the counter-jet direction. 
Possible explanations for the displacement include orbital motion of an SMBH binary, gravitational 
perturbations due to massive objects (e.g., globular clusters), acceleration by an asymmetric or in- 
trinsically one-sided jet, and gravitational recoil resulting from the coalescence of an SMBH binary. 
The displacement direction favors the latter two mechanisms. However, jet asymmetry is only viable, 
at the observed accretion rate, for a jet age of >0.1 Gyr and if the galaxy restoring force is negligible. 
This could be the case in the low density core of M87. A moderate recoil ~1 Myr ago might explain 
the disturbed nature of the nuclear gas disk, could be aligned with the jet axis, and can produce the 
observed offset. Alternatively, the displacement could be due to residual oscillations resulting from a 
large recoil that occurred in the aftermath of a major merger < 1 Gyr ago. 
Subject headings: black hole physics — galaxies: individual (M87) — galaxies: nuclei 



1. INTRODUCTION 

It is generally assumed that supermassive black holes 
(SMBHs) reside at the centers of their host galaxies. 
However, SMBHs can be significantly displaced from 
their central locations by asy mmetric forces during a 
merger, or by a second SMBH (|Komo ssa 2006). In addi- 
tion, if a binary SMBH forms and coalesces, anisotropic 
emission of gravitational waves can result in initial im- 
pulsive kick velocities of several thousand km s~^(e.g., 
iPretor ius 2007). Even if the kick velocity is small enough 
for the coalesced SMBH to remain in the galaxy, A^- 
body simulations have shown that the SMBH can os- 
cillate wit hin the bulge for ^ 1 Gyr before coming to 
rest (,Gualandris fc Merritt]l2008[ ). Alternatively, if the 
galaxy contains a radio source, the SMBH may experi- 
ence sustained acceleration due to intrinsic asymmetries 
in jet power (e.g., iTsygan 2007). 

The most direct way to find a displaced SMBH is to 
observe a spatial offset between the SMBH and the cen- 
ter of its host galaxy. This requires data at the high- 
est possible spatial resolution. E/SO galaxies that are 
minimally affected by extinction and contain an active 
galactic nucleus (AGN) provide good candidates; bulge 
isophotes determine the position of the galaxy center and 
the AGN (point source) determines the position of the 
SMBH. M87 is an ideal target for a displaced SMBH 
search. It is nearby, has a regular b ulge, is relatively fre e 
of dust, hosts an AGN and jet fe.g.. lPerlman et"ani2001[ ). 
and has been extensively observed by the Hubble Space 
Telescope [HST). In this Letter, we report the discovery 
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of a 6.8 ±0.8 pc projected displacement between the cen- 
ter of M87, as defined by the galaxy isophotes, and the 
SMBH. 

2. THE DATA 

Our analysis uses the archived HST data listed in Ta- 
ble [2 All data had the standard STScI on-the-fly re- 
processing applied. Each image was rotated to the same 
reference frame and shifted to a common position using 
a 2D cross correlation register. The data were combined 
using a median filter to remove residual cosmic rays and 
bad pixels, and to minimize the effect of the High Reso- 
lution Channel (HRC) coronagraphic aberration. 

The IRAF task ELLIPSE (Jcdrzejewski 1987) was used 
to determine the photo-center of the galaxy. Beginning 
with a semi-major axis (SMA) of 1 pixel, centered on the 
nuclear point source, ellipses of progressively increasing 
SMA were independently fitted to the data. The SMA 
was incremented by 1 pixel in each successive fit and the 
center of the ellipse found. The x — y pixel co-ordinates 
of the ellipse centers were thus determined as a function 
of SMA. To estimate the precision with which offsets can 
be recovered, and check whether masking degrades ac- 
curacy, we applied this technique to a set of simulated 
galaxies containing a nuclear point source. The simu- 
lated galaxies were given an r^''^ surface brightness pro- 
file and an ellipticity of 0.2. Point sources, with offsets 
of 0, 0.25, 0.5, 0.75, 1.0, 1.5, 2.0, 3.0 pixels in both x 
and y, were added to create 64 total models. Each model 
was convolved with the ffg r point sprea d function (PSF) 
as generated by TinyTim (|Kristi[l995l ). and populated 
with random noise. First, isophotal fits were performed 
without a mask. In this case, the offsets were recovered 
to ^0.2 pixels (Fig. [1^). Second, isophotes were fitted 
with a mask that simulated an extended jet crossing the 
galaxy center, multiple globular clusters and areas of ex- 
tinction. These offsets were also recovered to ~0.2 pixels 
(Fig. [1)3); any observed offsets will not be a result of 
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Fig. 1. — Accurate recovery of simulated SMBH-galaxy displace- 
ments, [a] Without [b] with a mask. Horizontal lines show the 
model displacements. Offsets are recovered for SMA<10 pixels. 
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Fig. 2.— la] ACS HRC F814W median combined image of M87 
with a logarithmic stretch, [b] The mask overlaid on [a]. 



TABLE 1 
Archived ACS Data 



TABLE 2 
Radial Offsets of Ellipse Centers 



Dataset 



Aperture 



Filter 



Exposure times 



PID 



J8q0* 
j92j* 
j9ei* 
j9qf* 
jSIO* 
jSqO* 
j92j* 
j9ei* 
j9qf* 
J9e0* 
J9e0* 



HRC 
HRC 
HRC 
HRC 
HRC 
HRC 
HRC 
HRC 
HRC 
WFC 
WFC 



F606W 
F606W 
F606W 
F606W 
F814W 
F814W 
F814W 
F814W 
F814W 
F606W 
F814W 



2 X 75s 
4 X 80s,5 X 48s 
6 X 80s,3 X 45s 
2 X 80s, Ix 45s 

5 X 100s 

8 X 50s 

4 X 50s 

6 X 48s 

2 X 48s 

3 X 500s 
3 X 1440s 



9829 
10133 
10617 
10910 
9705 
9829 
10133 
10617 
10910 
10543 
10543 



Note. — Details of the archived ACS data used. 
HST program ID number. 



PID is the 



masking. 

Figure [2] shows the median combined HRC data in 
F814W, both with and without the mask . A distance 
modulus of 31.0 ± 0.2 (iTonrv et all 1200 It ) puts M87 at 
16.1 ± 0.2 Mpc (77.9 pc per arcsec). The ACS pixel 
scales of 0'.'027 (HRC) and 0'.'049 (Wide Field Channel, 
WFC) are then 2.1 and 3.8 pc, respectively. 

3. RESULTS 

Table [2] presents the eUipse fits to the HRC and WFC 
data in F606W and F814W. The eUipticity and posi- 
tion angle (PA) of the fits are uncorrelated with SMA, 
show no evidence for b ulge asymmetry , and are consis- 
tent with the results of iFerrarese et al.l (|200(^ who note 
the isophotes are regular. At each SMA the radial off- 
set of the ellipse center is given with respect to the nu- 
clear point source. Figure [3] presents the x and y co- 
ordinates of the ellipse centers and the radial offset as a 
function of SMA. Note that the inner jet region is seen 
as "pear shaped" contours in Figure |3l^a,b). This feature 
was masked out during the fits. The elongated contours 
in Figure |3l(c,d) are an effect of charge bleeding from the 
strong point source, and result in larger uncertainties and 
the anonymously large offsets seen at ^I'.'O. 

All ellipses with SMA>1'.'0 show a clear offset relative 
to the nuclear point source. For SMA<1'.'0, the measured 
offsets show the transition between isophotes dominated 
by the AGN and the bulge. 

The mean radial offset from both the HRC and WFC 
(1'.'0<SMA<3'.'0), weighted by the uncertainties, is 6.84± 
0.07 pc. The error on the mean is substantially less than 
the 0.2 pixel uncertainty implied by the simulations (cor- 
responding to 0.76 pc in the WFC data), which we prefer 



HRC 
SMA (") Offset (pc) 

F606W F814W 



WFC 
SMA (") Offset (pc) 

F606W F814W 



0.027 


0.000 


0.000 


0.049 


0.000 


0.000 


0.054 


0.110 


0.187 


0.098 


0.122 


1.405 


0.081 


0.513 


0.453 


0.147 


0.602 


1.779 


0.108 


0.889 


1.155 


0.196 


1.272 


1.718 


0.135 


0.982 


1.350 


0.245 


2.071 


0.870 



Note. — SMA verses radial offset for all ACS data. Tabled is 
published in its entirety in the electronic edition of the ApJ Letters. 

to adopt as a conservative estimate of the offset uncer- 
tainty. Therefore, the best estimate of the radial offset is 
6.8 ± 0.8 pc. The uncertainty weighted PA between the 
offset and the p oint source is 306 i b 17°, consistent with 
the jet axis (e.g. lOwen et al.iri980[ ). 

4. ORIGIN OF THE DISPLACEMENT 

The observed offset between the bulge photo-center 
and the AGN suggests that the SMBH in M87 is dis- 
placed from the stellar center of mass. The projected 
displacement is ^ 7 pc (~ O'.'l) approximately in the 
counter-jet direction. Here we consider several mecha- 
nisms that might produce this displacement. 

4.1. Jet Asymmetry 

The PA of the disp lacement suggests a connection with 
the M87 jet. iShklov ski (1982) first noted that one-sided 
jets can accelerate SMBHs. However, the standard inter- 
pretation for one-sided jets is relativistic beaming (e.g., 
lEichler fc SmithI 119831 ). In the case of M87, there is 
VLBI and VLBA evidence for a p robable counter-jet 
()Lv et al.ll2007t iKovalev et"all 120071 ). the large-scale ra- 
dio structure shows evidence for two-sided jet activity 
in earlier epochs (Owen et al. 2000), and Chandra im- 
ages show symmetric X-ray inner cocoons (|Forman et al.l 
l2007f ) that could be used to constrain the jet asymmetry. 
Models in which the SMBH is accelerated by two-sided, 
but intrinsically as ym metric, jet s have been explored by 
iWang et al.l (|1992D and lTsveanl ([2007). In this case, the 
SMBH acceleration is 



obh 



: 2.1 X lO'Vi 



fjctTO cm s ^ (1) 

(jKornreich fc Lovelacel[2008l) where /jet is the luminosity 
of the asymmetric part of the jet, expressed as a fraction 
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Fig. 3. — [Top row] x and y ellipse centers overlaid on ACS contours. The point source is at [0,0]. [a,b] HRC ellipse centers for SMA<1'.'5. 
The jet (masked out for the fits) can be seen as a N-W extension. [c,d] WFC ellipse centers for SMA<3'.'0. {a,cj F606W {h,dj F814W. 
[Bottom row] Radial offsets of the ellipse centers. 



of the accretion luminosity. The mass accretion rate (m) 
is in units of AfEdd = LEdd/{^c^), where e is the accretion 
efficiency. 

As M87 has a large, low- density core (e.g.. lLauer et al.l 
Il992t iFerrarese et al.ll2006f ) , we first assume the restoring 
force from the galaxy is negligible. The displacement 
(Ar) and velocity (Av) of the SMBH would increase with 
time as 



In the presence of a restoring force, the SMBH is ex- 
pected to come to rest where the force from the galaxy 
matches the jet force. Adopting units such that G — 
Tc = (T — 1, where re ~ 500 pc is the M87 core radius 
and a « 330km s^^ is the Id core stellar velocity disper- 
sion, the equation of moti on of an accelerated SMBH in 
a fixed galaxy core is (e.g. iGualandris fc Merrittll2008l ) 



Ar«340pc/jot™t^, 
Av « 660km s^ /jctw ^6 



(2a) 

(2b) 

where tg is the time in Myr since the jet turned on. As- 
suming the SMBH is offset in the counter-jet direction, a 
projected displa cement of 6.8 ib 0.8 pc, and a jet orienta- 
tion of 15 ± 5° ([Biretta et al .111 9991 ). gives a physical Ar 
of 26lg® pc. However, this could be as small as ^ 10 pc 
for a jet orientation of 45° (|Lv et al.ll2007D . An upper 
limit to the jet li fetime is set by ag e of the outer radio 
halo {^ O.lGvr: lOwen et aLll2000D . assuming it is stiU 
powered by the jet. Alternatively, if the outer halo is 
a relic, then a lower limit to the jet li fetime is given by 
the a ge of the inner lobes (~1 Mvr: iBicknell fc BegelmanI 
119961) . Equations ^ therefore give 

3xl0-6(i^)</jetm<3xl0-2(^^), (3a) 
0-2 {^) km s-i < At; < 20 (^) km s-^. (3b) 



Xi + 2jXi + [3'^Xi = Qi 
where the cc, are Cartesian coordinates and 
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(5a) 

(5b) 
(5c) 
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iDi Matteo et al.l (2003) report m « lO"'' in M87. A 
1 Myr jet then requires /jet >> 1 and can be ruled 
out. However, for a ~0.1 Gyr jet, fjet ~ 0.03, i.e., the 
jet asymmetry amounts to only ~ 3% of the accretion lu- 
minosity. Therefore, the displacement could result from 
acceleration by a long-lived jet if the galaxy restoring 
force is negligible. 



The third term on the LHS of equation 2] represents the 
restoring force on the displaced SMBH from the stars, 
assuming that the motion takes place in a homogeneous 
core; the Ci are related to the sh ape of the core and ar e 
unity for a spherical galaxy (e.g.. lChandrasekhaHll987D . 
The second term represents the dynamical friction force 
from the stars and has been expressed in terms of the core 
mass, Mc = Anprc^ /3; In A is the Coulomb logarithm and 
F <; 1 is a "fudge factor" accounting for the fact that the 
frictional force on massive objects in galaxy cores is found 
to be soi newhat less than predicted by Chandrase khar's 
formula (jGualandris fc MerritH [200l Ihlmii l200l . Fi- 
nally, Xi = Cj ■ Si where Cj is a unit vector in the direction 
of the jet, and r5oo = rc/(500pc), 0-300 = cr/(300km s"-^). 
In the absence of dynamical friction (7 = 0), the solu- 
tions to equation ^ are 



Xi{t) = —Xi {cosujit — 1) , 



X,, 



4, (6) 



w; 



Batcheldor et al. 



i.e. continued oscillation about the point X where the 
jet acceleration is balanced by the gravitational force. 
Unless the core is very non-spherical, X will point ap- 
proximately opposite the jet direction. 

In the presence of dynamical friction, the oscillations 
are damped and asymptote to X. Furthermore, since 
M, is of order Ale, the damping time is expected to 
be comparable to the core crossing time, Tc = rc/<Jc ~ 
1.5 Myr {r^oo / <^30o) ■ This is somewhat larger if F < 1, 
but is nevertheless comparable with the lower limit on 
the jet lifetime set by the age of the inner lobes. 

Reintroducing dimensional variables, 

Ar ~ 500pc ifjctm) r^oQf^aoo- (7) 

This implies 0.02 <, fj^tm <, 0.05 for 
lOpc ^ Ar ^ 25pc. Even assuming /jet — 1 (one- 
sided jet), the observed offset requires m sa 0.02, i.e., 
several magnitudes more than observed. In addition, 
while the AGN in M87 may have been more active in 
the past, the offset requires the jet to have remained at 
this high level of activity for >1 Myr. 

4.2. Binary SMBH 

If the SMBH in M87 is one component of a bound pair, 
with mass ratio q = AI2/M1 < 1, then 



motion of the SMBH, with mean square velocity 



qAr2 



Mil) 



Vioopcy 



(8) 



where Ar2 is the separation of the smaller SMBH from 
the binary center of mass (located at the galaxy center). 
A binary separation of ~ 100 pc is not unreasonable given 
the high merger rate expected for a luminous galaxy at 
the center of a rich cluster, and given the gradual rate of 
in-spiral expected in the low-density core. 

The velocity of the larger SMBH with respect to the 
binary center of mass (assuming a circular orbit) would 
be much greater than in the jet scenario: 
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G(Afi+M2) 



400km 
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1+9 



f M1+M2 y^^ (" g A 
\4xl0'-> Mq J l^llOpc^ 



-1/2 



(9a) 
(9b) 



where a is the binary separation. However, the line-of- 
sight velocity is Vjos = Vi sin i sin 0, where i and (p are 
the unknown inclination and phase of the orbit. If the 
nuclear gas disk is in the orbital plane then i w 50° 
(jMacchetto et all I1997D which gives Vios ~ 220km s~i 
when (j) = 45°. 

A non-active second SMBH would be extremely hard 
to detect, but a likely consequen ce of such a binary would 
be jet precession (e.g., Romero et al.ll2000h . Several jet 
knots have strong helical morphologies, but their widths 
are consistent with a steadi ly expanding cone wi th near- 
zero width at the nucleus (jLobanov et al.|[2005[ ). Since 
there is no evidence of a O'.'l scale "wobble" of the jet 
direction in both HST and VLBA data, any precession 
must be several orders of magnitude smaller than that 
expected from a binary. 

4.3. Massive Perturbers 

Any SMBH will experience gravitational perturbations 
from stars and more massive objects, e.g. globular clus- 
ters (GCs), open clusters, etc. The result is a Brownian 



-Af,(VsMBH) 



-ma . 



(10) 



Here, fh is the second moment of the mass distribution of 
perturbing objects, and a is the velocity dispersion of the 
perturbers meas ured within ~ O.Sr ,; (r^ is the SMBH's 
influence radius. iMerritt et al.ll2007l ). The rms displace- 
ment of the SMBH is then given by the virial theorem: 



Ari-ms w — - 



m 



1/2 



(11) 



Assumin g that GCs constitu te a fraction ^ 10^"^ by mass 
of M87 iMcLaughlin et all [1994) and that the mass of 
one GC is lO^M©, Equation [TT] implies rims ^0.1 pc. 
This is likely an overestimate since the density profile 
of GCs is flatter than that of the galaxy light. Unless 
there is another popula tion of "massive perturbers" in 
M87 (jPerets et al.ll2007[ ) it is unlikely that the observed 
displacement can be due to gravitational perturbations. 

4.4. Gravitational Wave Recoil 

A kick velocity of Uk ~ 500km s~^ would displace the 
coalesced SMBH a distance ~ re from the center of M87, 
in a direction orthogonal to the orbital plane of the pre- 
ceding binary (|van Meter et al.|[2blQ ). and aligned with 
the jet axis. Such a velocity is easily produced during 
the coalescence of t wo modestly spinning, com parably 
massive SMBHs (e.g.. lTichv fc Marronettill2007[ ). Larger 
kicks (but smaller than the escape velocity) would result 
in SMBH-core oscillations that damp on a time scale 



T, 



damp 



15 



G2pAf. 



2 X 10«yr ^aor'^soo (12) 



(jGualandris fc Merrittl [20081 ) . The rms displacement of 
the SMBH with respect to the galaxy center decreases 
with time as 



-m 



^(t-te)/2Td„ 



t > tr. 



(13) 



where tc is the time at which the amplitude of the 
oscillations has damped to a scale of '^ Tc. A cur- 
rent displacement of ~ 20 pc ~ 0.05rc implies that 
^ STdamp ~ 10 Gyr have elapsed since tc- Thus, a major 
SMBH merger during the formation of M87 could have 
generated a kick consistent with the current displace- 
ment. 

5. DISCUSSION 

Four displacement scenarios have been considered; jets, 
binaries, perturbers, and gravitational recoil. Perturbers 
cannot produce the observed displacement amplitude, 
and the observed jet-offset alignment must be fortuitous 
in the binary case. On the other hand, jet acceleration 
and gravitational recoil provide natural explanations for 
the jet-offset alignment. 

These mechanisms predict displacement velocities that 
are either <^ 10km s~^ or ^ 100km s~^. For example, 
jet acceleration of the SMBH produces an rms velocity 
Wrms ~ 0km s~^ in the presence of a restoring force. 



and Aw 



20km s~^ t7^ without. For SMBH-core 
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oscillations induced by a large and early kick, the rms 
velocity of the SMBH is 



Wrr 



r,.,; 



-101 -1 —l^rms 

12kms 0-300 rgoo: 



'20 pc' 



(14) 



independent of the time since the kick or the kick am- 
plitude. However, a velocity of several hundred km s~^ 
is possible if the SMBH is undergoing a large amplitude 
oscillation following a recent (< 1 Myr) kick. A similar 
velocity is expected in the binary SMBH scenario. 

It may be possible to distinguish between the high 
and low velocity displacement mechanisms by comparing 
the relative recessional velocities of stars at the galaxy 
center and gas around the SMBH. Unfortunately, this 
will be difficult; the low central surface brightness will 
limit the accuracy with which velocities can be deter- 
mined from stellar absorption lines. Alternatively, wide 
field VL BA astrometry, ref erenced against background 
quasars (jDavies et al.ll2009() . may allow the proper mo- 
tion of the nuclear point source to be determined. How- 
ever, it will be difficult to distinguish between the cluster 
motion of M87 and the intrinsic kinematics of the SMBH. 

Th e central gas disk (Harms et al.l 119941 : [Ford et aH 
II 9941) may provide alternative constraints on the dis- 
placement velocity. For example, a displaced SMBH re- 
tains g as whose Keplerian v elocity exceeds th e kick ve- 
locity (jMerritt et al.ll2006D . iFord et al.l ()1994[ ) estimate 
the disk-like structure in M87 to extend to ^ 80 pc. 
Assuming this represents the portion of the disk that 
remains bound to the SMBH, we can infer an initial 
kick velocity '^ 400km s~^. This would produce strongly 
shocked regions in the disk, as it passes through the 
ambient gas, that may explain several observed proper- 
ties. For example, while the gas in the central ifc0'.^4 re- 



gion of the disk is p hoto-ionized (jMacchetto et al.lll997l : 
iSankrit et al.lll99^ and has a closely Kepleria n velocity 
field (jHarms et al.lll994tlMacchetto et al.lll997| ). at larger 
radii non circul ar motions and morphological distor- 
tions are present (jFord fc Tsvetanovll 19991 iBradlev et al.l 
[2004) . In addition, gas in the outer part s of the disk 
is evi dently excited by a 265km s^^shock (jDopita et al.l 
Il997f) . rather than the central UV continuum source. 

Two further consequences from a kick velocity t^k ~ 
400km s~^ are noted. First, the shocked gas may be 
hot enough (T ^ 2 x 10^ K) to produce X-ray emis- 
sion. Second, at this velocity the observed displacement 
is reached after ~ 25000(Ar/10pc) yr, i.e., significantly 
less than the minimum jet age discussed earlier. If this 
minimum jet age corresponds to the time since the kick. 



the SMBH cannot be on its initial outward trajectory. 
iGualandris fc MerrittI (|2008[ ) show for Vk < 50% of the 
escape velocity, the SMBH motion is rapidly damped and 
undergoes only one large amplitude oscillation on a time- 
scale of 1 Myr . Therefore, it is possible the SMBH is on 
a return trajectory. 

6. SUMMARY 

We find that the SMBH in M87 is displaced relative 
to the galaxy photo-center by a projected distance of 
6.8 ± 0.8 pc (^ O'.'l) in the counter-jet direction. Four 
explanations have been considered: jet-induced acceler- 
ation, a binary SMBH, massive perturbers and a grav- 
itational wave kick. Perturbers only produce offsets of 
^0.1 pc, and precession of an SMBH binary would pro- 
duce "wiggles" in the jet that are not seen. Neither can 
explain the observed alignment with the jet axis. Jet 
acceleration cannot be ruled out, but it requires both 
that the jet age to be >> 1 Myr and that the restoring 
force exerted by the galaxy is small. The displacement 
could also be be explained by a moderate (a few hundred 
km s~^) kick which occurred ^ 1 Myr ago. This could 
produce the alignment with the jet axis, and may also 
explain the disturbed nature of the nuclear gas disk. Al- 
ternatively, the displacement may be from SMBH-core 
oscillations following a kick that occurred <1 Gyr ago. 

The displacement processes discussed are likely to be 
common in early type galaxies; they are partly assembled 
via mergers and often host AGN-powered radio sources. 
It is plausible, therefore, that there is a high incidence 
of SMBH displacements among the E/SO population. A 
systematic effort to determine the statistical distribution 
of displaced nuclei will provide important insights into 
merger histories, the frequency with which SMBH bina- 
ries form and coalesce, and/or jet launch physics. What- 
ever the mechanism, displaced SMBHs have important 
consequences for accretion flows, for stellar dynamics in 
galactic nuclei, and for our fundamental understanding 
of galaxy evolution. 
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